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ABSTRACT 

Multi-wavelength imaging polarimetry at far-infrared wavelengths has proven to be an excellent tool for studying 
the physical properties of dust, molecular clouds, and magnetic fields in the interstellar medium. Although these 
wavelengths are only observable from airborne or space-based platforms, no first-generation instrument for the 
Stratospheric Observatory for Infrared Astronomy (SOFIA) is presently designed with polarimetric capabilities. 
We study several options for upgrading the High-resolution Airborne Wideband Camera (HAWC) to a sensitive 
FIR polarimcter. HAWC is a 12 x 32 pixel bolometer camera designed to cover the 53 - 215 fim spectral range 
in 4 colors, all at diffraction-limited resolution (5 - 21 arcsec). Upgrade options include: (1) an external set of 
optics which modulates the polarization state of the incoming radiation before entering the cryostat window; 
(2) internal polarizing optics; and (3) a replacement of the current detector array with two state-of-the-art 
superconducting bolometer arrays, an upgrade of the HAWC camera as well as polarimeter. We discuss a range 
of science studies which will be possible with these upgrades including magnetic fields in star-forming regions 
and galaxies and the wavelength-dependence of polarization. 

Keywords: polarimetry, far-infrared astronomy, instrumentation, interstellar medium, magnetic fields, dust, 
airborne astronomy, SOFIA 

1. INTRODUCTION 

Mapping polarimetry at far-infrarcd/submillimeter wavelengths began on the Kuiper Airborne Observatory 
(KAO) more than 20 years ago. Over the course of the KAO's lifetime polarimetric instruments evolved from 
single pixel devices 1 to larger arrays 2 and eventually moved to ground-based observatories. 3,4 Far-infrared (FIR) 
and submillimeter observations of dense molecular clouds, protostellar disks and envelopes, the Galactic center, 
and one external galaxy have shown that dust emission from these objects is polarized (at levels of ~ 2 - 10%) at 
almost every point. 5 ~ 7 The alignment of dust grains with the local magnetic field 8 allows one to infer the direction 
of the magnetic field projected onto the plane of the sky. This allows studies of magnetic fields' interactions with 
the local interstellar medium and their role in star and galaxy formation 9-12 as well as studies of dust and cloud 
properties. 13, 14 

As successor to the KAO, the Stratospheric Observatory for Infrared Astronomy (SOFIA) offers a new 
frontier for observations in the infrared-to-submillimeter spectral regime. Unfortunately, no first-light SOFIA 
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Figure 1. The HAWC cryostat showing the locations of several proposed polarimetric upgrades: dual-polarization detector 
arrays, wave-plates (HWP), and the variable-delay modulator. 



instruments 15 have existing capability for accurate polarimetry. Both a far-infrared (Hale) 16 and mid-infrared 17 
polarimeter have been proposed as future SOFIA instruments. Here we discuss three options for upgrading the 
High-resolution Airborne Wideband Camera (HAWC) 18 to a sensitive FIR polarimeter. 

HAWC is a first-light instrument for SOFIA. It consists of a 12 x 32 pixel bolometer camera designed to cover 
the 53 - 215 /mi spectral range in 4 colors, all at diffraction- limited resolution (5" - 21"). A complete upgrade 
of the HAWC camera into a polarimeter (which we call Hale) will include replacing the existing 384 pixel array 
of silicon bolometers with two superconducting bolometer arrays. Each such array would consist of over 5000 
pixels and detect orthogonal polarization components selected by wire grids close to the detectors. A cold (4K) 
crystal half-wave plate (HWP) upstream from the grids modulates the incoming radiation. In this way we create 
a sensitive, diffraction- limited, 4-color polarimeter. Additionally, the dual-polarization measurement and the 
detector upgrade yield a photometer with a much larger field-of-view than the original HAWC instrument. 

We consider two other options (called HAWC-pol) as "pathfinder" programs working towards the more 
complete Hale instrument. These options are attractive as it is feasible that they may be incorporated into 
HAWC soon after its first light on SOFIA. In both HAWC-pol options we retain the existing 12 x 32 silicon 
bolometer array and measure only a single component of polarization. Option (1) places one or more rotating 
HWPs at the cold HAWC pupil stop, followed by a polarizing grid. In option (2) we place only a polarizing grid 
at the pupil and modulate the polarization outside the cryostat using a variable-delay polarization modulator, or 
polarization "switch". Details of the upgrade components are discussed in Sect. [3] their locations in the HAWC 
cryostat are shown in Fig.[l] The technical specifications for both Hale and HAWC-pol are shown in Table [l] 



J. E. Vaillancourt et al. 



3 





Table 1. HAWC-pol and Hale 


specifications. 
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"Noise Equivalent Power: background limited, per pixel 

''Noise Equivalent Flux Density for HAWC camera: background limited, per beam, chopped 

c In 5 hours for 5 Jy source; chopped; assumes 100% polarization efficiency and 100% observing efficiency. 



2. SCIENTIFIC GOALS 

2.1. Turbulent Star Formation 

The paradigm for star formation, once seen as depending on slow diffusion of magnetic fields out of cloud cores, 
has shifted to a violent one, 19 the dynamics of which are governed by magnetic compressible turbulence. 20,21 

A unique view of the turbulent field is provided by polarimetry of emission from magnetically aligned grains. 
Problems of interpreting the statistics of 2D polarization maps in terms of the underlying statistics of the 3D 
magnetic fields are becoming tractable with modern theoretical work on radiative transfer, grain alignment, and 
turbulence, even in the case of inhomogeneous clouds containing dense clumps and clusters of embedded stars. 22 
The problem of field bending due to large-scale non-turbulent effects is being addressed by choosing reference 
systems along the local projected direction of the field rather than the mean field for the whole cloud. 

When maps of polarized emission can be made with SOFIA'S resolution, dynamic range, and spectral coverage, 
it will be feasible, not only to test predictions of turbulent vs. static star formation, but also to determine 
characteristics of the magnetic fields such as the ratio of fluctuating to mean components, the energy spectrum, 
and the distribution of magnetic field strengths. 

2.2. Dust Grains and the Polarization Spectrum 

Recent years have been marked by substantial progress in the understanding of grain alignment. 23 Research on 
radiative torques 24-27 and grain dynamics 28-35 have made it possible to explain the otherwise perplexing results 
of FIR polarization from starless cores 36 and to make quantitative predictions of grain alignment along lines of 
sight through dense clouds. 

The anisotropic radiation responsible for aligning grains is ubiquitous in astrophysics. Thus polarized emission 
from dust grains can trace magnetic fields in various astrophysical environments. (While radiative torques 
provide the aligning mechanism, the alignment occurs with respect to the local magnetic field.) Features in the 
measured polarization spectra of molecular clouds (Fig. [2]) have been attributed to varying dust temperatures 
and polarizing efficiencies along the line-of-sight. 13 ' 37 ' 38 The alignment of grains, their temperature, and their 
emission properties depend on the radiation field to which they are exposed. As a result, multi-wavelength 



polarimetry can be used for detailed studies of magnetic fields even within unresolved objects (e.g. Sect. 2.3 
Fig. [2). 

A polarimeter on SOFIA will make it possible to pursue these studies with more passbands, more pixels, and 
improved spatial resolution. 
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Figure 2. Left: Measured polarization spectra of several molecular cloud envelopes, each cloud is normalized at 
350/im. 38 ' 39 Right: Predicted spectral energy distribution (bottom) and polarization spectrum (top) for protoplane- 
tary disk emission. 40 Note that the peak in the polarization occurs at 45 /jm, corresponding approximately to the shortest 
HAWC passband. 



2.3. Polarization spectra of protoplanetary disks 

By the time a young stellar object reaches the Class II stage, a.k.a. T-Tauri star (TTS) stage, the star has acquired 
approximately its final mass, but a residual disk still contains diffuse matter that will presumably be incorporated 
into planets. According to some current theories of planet formation, 41 the first step is for the grains in such a 
protoplanetary disk to coagulate, forming particles with sizes of order one millimeter or larger. Observational 
evidence of such "large grains" has been obtained via modeling of the submillimeter-millimeter spectral energy 
distributions of TTS disks. 42-45 If these models are correct, then the optical depths for absorption and scattering 
will approach or exceed unity across much of the FIR-to-millimeter spectrum. 43 In this case, the polarization of 
the disk thermal emission will likely contain contributions from selective absorption and scattering. 46 Although 
this complicates the interpretation, it also presents us with a new opportunity: Because the cross-sections for 
absorption and scattering, and thus the polarization characteristics, depend strongly on grain size for sizes 
comparable to the wavelength, we may be able to confirm the existence of millimeter-sized grains and further 
constrain the grain size via polarimetry of the disk emission. 

Polarimetry of protoplanetary disk emission has been sparse and the modeling is at an early stage, but with 
the advent of SOFIA, ALMA0 and EVLA 47 it will become possible to carry out sensitive polarimetry across the 
FIR-to-millimeter wavebands for a sample of bright TTS disks. A recent model 40 shows that the polarization 
characteristics depend strongly on grain size (Fig. |2j. Initial polarimetric observations of TTS disk emission, 
obtained at the James Clerk Maxwell Telescope 48 and at the Caltech Submillimeter Observatory (M. Krejny, 
priv. comm.) suggest polarization levels of one to several percent. Four-color polarimetry with Hale or HAWC- 
pol will place constraints on grain size and provide information on the magnetic fields (the latter of which may 
influence disk evolution 49 ). 



* http: / / www.alma. nrao.edu/ 
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2.4. Connecting Large- and Small-Scale Interstellar Fields 

The role played by magnetic fields in star formation has remained uncertain. At one extreme of the theories of 
star formation is that magnetic fields control the formation and evolution of the molecular clouds from which stars 
form, including the formation of cores and their gravitational collapse to form protostars. 50 The other extreme is 
that magnetic fields are unimportant, with molecular clouds forming at the intersection of turbulent supersonic 
flows in the interstellar medium. 51 In spite of considerable effort to observe magnetic fields in molecular clouds in 
order to resolve this uncertainty, the issue remains controversial. SOFIA with Hale or HAWC-pol will provide a 
crucial link between large- and small-scale data on interstellar magnetic fields that will be essential to understand 
the role of magnetic fields in star formation. 

On the largest scales, observations of stellar polarization 52 ' 53 and Planck 54 polarization maps at 5' resolution 
enable mapping of the Galactic to molecular cloud scale structure of the interstellar magnetic field. HAWC-pol 
will make it possible to study structure of cores relative to large-scale magnetic fields. Some predicted phenomena 
that require the correlation of the large- and small-scale maps of magnetic field structures include collapse of 
mass along field lines to form cores flattened along field lines, hourglass morphology fields in cores with the core 
fields connecting to the larger scale molecular cloud fields, magnetic braking of cores that will twist the fields as 
angular momentum is transferred outward from cores to envelopes, and bipolar outflows from protostars with 
magnetic fields being parallel to the outflows. 

On the smallest scales (< few arcseconds), ALMA will image protostars, protoplanetary disks, and the 
inner parts of protostellar outflows. In spite of the much larger telescope collecting area of ALMA, the shorter 
wavelengths and much higher bandwidths observed by HAWC mean that HAWC-pol will have comparable 
sensitivity to ALMA for polarization mapping of extended emission. Theoretical models predict the morphology 
of the connection of magnetic fields in these very small scale phenomena to the larger scale fields in the cores 
from which protostars form. Combining ALMA and HAWC-pol maps will probe the full spatial structure of 
these phenomena and allow testing of the star formation theoretical results. 

2.5. The Galactic Center and External Galaxies 

Our proximity to the Galactic Center provides an opportunity to study the physics of galactic nuclei in great 
detail. There is evidence that magnetic fields play an important role in the dynamics of the Galactic center region 
as evidenced by the prominent Radio Arc 55 and other non-thermal filaments 56,57 that exist in this region. These 
filaments are thought to be the result of relativistic electrons spiraling along lines of magnetic flux, and though 
the brightest arc oriented perpendicular to the plane, more recently discovered filaments have been found with 
other orientations. Submillimctcr polarimetry has indicated that the large-scale magnetic field in the molecular 
gas is generally parallel to the plane 58 in stark contrast to the brightest radio filaments. However, on smaller 
scales the field has been found to be more complex (Fig. [3}. 10 

Key outstanding questions concerning the magnetic fields in the Galactic center are: (1) What is the geometry 
of the field, (2) what is the strength of the field, and (3) how do the radio filaments form? The multi-frequency 
capability of a HAWC polarimeter will allow us to separately measure the field in different components along 
the line-of-sight. This is particularly important in the Galactic center, as the region encompasses a wide range 
of temperatures. In addition, the high angular resolution may enable one to probe the interaction sites between 
filaments and associated molecular clouds. Such interactions may provide a mechanism for electron acceleration 
which is required for the formation of the filaments. This has the potential for providing estimates of the local 
magnetic field strength and constraining models for filament formation. 

External galaxies give us the opportunity to study the magnetic field geometry of an entire system viewed 
from the outside. Edge on spiral galaxies in particular will allow us to sample the integrated interstellar medium 
along a number of lines of sight, weighted of course by the distribution of dust, stars, and heating radiation as a 
function of galactocentric radius. For many galaxies such as NGC891, NGC4565, etc., the edge-on disk is only 
about one beam wide at 60 /jm and spans many beam widths across the sky. Comparison of very active galaxies 
such as M82, moderately active disks such as in NGC891 and more quiescent disks such as in NGC4565 will 
allow Hale to find and distinguish large-scale structures in the magnetic field deep within these galaxies. These 
features could include massive blowouts, for example, with a projected magnetic field geometry vertical to the 
galactic plane, undetectable with optical and near-infrared polarimetry. 
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Figure 3. Polarization and radio map of the Galactic center. 10 Inferred magnetic field vectors from FIR/submillimeter 
polarimetry are superposed on a 20 cm continuum VLA image. The 100 fim vectors appear to trace the thermal arched 
filaments while the magnetic field in the molecular cloud associated with GO. 18-0.04 (350 fim polarimetry) is perpendicular 
to the field traced by the non-thermal filaments of the Radio Arc. 



3. INSTRUMENT DESIGN 

3.1. Polarization Modulation 

Atmospheric constraints on suborbital FIR/submillimeter polarimetry have been described in the literature. 59,60 
The two main effects are: 1) the temporal change of the atmospheric transmission, and 2) "sky noise" , which 
is the spatially and temporally variable emission from the atmosphere which is not completely eliminated by 
techniques such as chopping. Both effects can be mitigated by a polarimeter which observes the field of view 
in two polarizations simultaneously. 60 However, it is often not practical to implement the dual-polarization 
approach in an existing instrument design. An alternate solution is to perform polarization modulation quickly, 
so as to "freeze" the atmosphere. 59 ' 61 Both the dual-polarization and rapid-polarization-modulation techniques 
work because atmospheric emission is intrinsically unpolarized. 

For HAWC-pol, we are considering two approaches to rapid polarization modulation - a variable-delay po- 
larization modulator (VPM), and a continuously-spinning half-wave plate (HWP). Both devices are considered 
in more detail in the following sections. However, several aspects of the design are common. 

In our observing plan, the chopping secondary mirror of SOFIA is the fastest modulation, and the VPM or 
HWP is a slower modulation. Secondary chopping is ideally done at a suitably high frequency where the noise 
of the difference is dominated by white noise from the photon arrival statistics. 62 The typical atmospheric noise 
spectrum for SOFIA FIR observations is not yet known, but the mechanical capabilities support chopping at 5 
- 20 Hz, which we expect will be fast enough in most conditions. 

Although the atmospheric sky noise is intrinsically unpolarized, it will become partially polarized by non- 
normal reflections in the optics, in particular the dichroic tertiary mirror. Past experience with this type of 
mirror 5 indicates that dichroic mirror polarization will be ~ 1%. Therefore, the sky noise will be registered as 
a polarization noise at a level ~ 100 x reduced. Since the polarization modulation is at a lower frequency, the 
frequency spectrum of the atmospheric noise matters. If the noise spectrum is proportional to l// 2 (in amplitude 
units), then the polarization modulation needs to be at a frequency at least as large as ~ O.lx the requirement 
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Figure 4. Left: A schematic for the Variable-delay Polarization Modulator (VPM). A controlled, variable delay (shaded 
line) can be introduced between two orthogonal linear polarizations. Right: The current HAWC fore-optics can be 
modified for the polarimeter. The optics to the left of the window are at the ambient temperature; these include the 
VPM and the Focusing Mirror (FM). The optics to the right of the window are all ~ 4K; these include the analyzing 
grid, filters, lens, detector folding mirror (DFM), and the detectors themselves. 



for the chopping. For HAWC-pol, we plan on polarization modulation at ~ 1 Hz with either a square- wave profile 
(VPM) or sinusoidal profile (HWP or VPM). 

Since all of the radiation from the atmosphere is partially polarized by the dichroic mirror, continuous 
polarization modulation will produce a large, nearly constant periodic signal on the detectors equivalent to a 
~ IK (Rayleigh- Jeans) load. The detector readout must have a dynamic range which accepts this signal, and 
the linearity of the detector over this range must be understood. Since the polarization modulation is used with 
faster spatial modulation, the celestial polarization signal will appear at ~ 1Hz "sidebands" of the chopping 
frequency and its harmonics; the detector readout must have good noise performance in this frequency range. 
All of these requirements are met by the HAWC system. The large 1 Hz atmospheric signal can be filtered out 
in an early data analysis step. 

3.2. A Continuously-Spinning Half- Wave Plate 

In this section, we describe an approach to HAWC polarimetry using a uniformly rotating bircfringcnt quartz 
HWP. Cryogenic crystal quartz has well characterized o-ray and e-ray indices of refraction, and thin wave- 
plates should have loss < 5% over most of the FIR. 63,64 Although waveplates can be stacked to achieve near- 
achromaticity over an octave of wavelength, 4 the HAWC filter bands cover two octaves. We expect that at least 
two HWPs will be required and perhaps four. 

We have identified the HAWC pupil wheel as an ideal location to locate polarimeter HWPs. An existing 
mechanism provides infrequent selection of up to eight pupil apertures for diagnostic purposes. In practice, no 
more than four of these apertures are needed for the camera, 18 leaving four positions available for HWPs. The 
pupil wheel is located at an image of the SOFIA primary mirror, which is an optimum location for a HWP. 65 
The pupil image has a diameter of 38 mm, a manageable size for quartz elements. A wire-grid polarizer will be 
located behind each HWP to provide polarization detection. The polarimeter is removed from the optical path 
by rotating the pupil wheel to select the camera-mode pupil aperture. 

For 1 Hz polarization modulation, the HWP must rotate at 15 rpm. Even at this low speed, induced vibration 
and microphonic response of the high-impedance bolometer detectors is a concern. We will avoid ball bearings 
and gear trains in the design and will instead use jewel bearings and three rollers to suspend each half-wave plate 
and a motor which directly torques permanent magnets attached to the waveplate. We expect dissipation in the 
bearings to be under 1 mW, which can be carried away through the pupil wheel. An optical encoder will monitor 
the HWP position. Implementing this polarimeter solution is straightforward and will require reconstruction of 
the pupil wheel, the plate which holds it, and the Geneva mechanism which defines the eight positions. 

3.3. The Variable-delay Polarization Modulator (VPM) 

In the case of the Variable-delay Polarization Modulator (VPM) , polarization is modulated by allowing a con- 
trolled phase delay to be introduced between two orthogonal linear polarizations. The basic principle is shown 
in Fig. HI and consists of a grid/mirror pair. Incident partially-polarized radiation comes in from the left. One 
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Table 2. Measuring Polarization with a VPM 
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component of polarization is reflected from the front of a grid. The other passes through the grid to a parallel 
mirror and is reflected back through the grid where the two beams are recombined. The component transmitted 
by the grid experiences an additional path length (which is shaded in the figure). As the grid-mirror separation is 
modulated, the polarization measured by polarization-sensitive detectors changes in a predictable way. 66 Various 
instruments have utilized devices with this architecture and have referred to them by various names. 66-70 The 
nomenclature "VPM" used here emphasizes the functionality as a polarization modulator and that the degree 
of freedom employed is that of a variable phase delay between linear polarizations. 

The polarization modulation for HAWC-pol must be flexible enough to accommodate all 4 HAWC passbands. 
As long as the grid is sufficiently fine for the 53 /xm band, we only need to change the grid-mirror separations to 
optimize the VPM for polarization modulation across the different bands. 

By making two changes to the current HAWC configuration, we can use a single VPM to measure Stokes Q 
and U for the four HAWC passbands: 1) Replace the warm folding flat mirror with a VPM that is mounted on a 
rotatable stage, and 2) Add two analyzer grids to the open slots on the aperture (pupil) wheel that have a relative 
orientation of 45°. Figure [4] also shows the warm and cold optics for HAWC with the proposed modifications for 
a VPM-based polarimeter. 

The polarimeter will have two modes, one to measure Stokes Q in the instrument coordinate system and 
the other to measure Stokes U. These will be differentiated by use of the two different analyzer grids that are 
oriented at a 45° with respect to one another. The VPM is placed on a rotator such that the relative angle of 
the VPM wires with respect to whichever analyzer grid is used is always 45°. In this way the measured signal 
on the detectors in each of the modes is 

S Q = I + Qcos(A^) + ^sin(A0), (1) 
SV = /+^cos(A0) + Usin(A0), (2) 

where A(f> is the phase delay introduced by the VPM. This can be scaled to any HAWC frequency by changing 
the amplitude of the grid-mirror separation. By using half- and full-wave delays for a given band, we can use 
the VPM as a polarization "switch" to rapidly change the polarization state to which the detectors are sensitive 
from I+Q — ► I — Q or I + U — > I — U depending upon the mode of operation. Table [2] illustrates the polarization 
measurement strategy. 

3.4. New Technology Large-format Detector Arrays 

To enable revolutionary FIR polarimetry, Hale will require simultaneous detection of two orthogonal linear 
polarizations. To accomplish this, we plan to implement two focal planes with 5,120 detectors each, using 
Goddard Space Flight Center's Backshort-Under-Grid (BUG) technology. 71 As in the SCUBA-2 design, 72 four 
such 32 x 40 subarrays will be tiled to produce each 5,120 element array. 

The BUG architecture is a two-piece, keyed assembly in which a two-dimensional array of backshorts are 
mounted under a matching array of suspended membranes. (It should be noted here that for Hale, the term 
"backshort" is a misnomer because in order to enable good performance over all four passbands, we will actually 
employ a "back termination."). Transition-Edge Sensors (TES) are fabricated on the membranes and the leads 
are "wrapped around" the frame to the back of the device at which point electrical connections to the multiplexer 
can be made via bump-bonding. The pixel size in this architecture is 1.135 mm so that an array of BUGs can 
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Figure 5. a) The Backshort-Under-Grid (BUGs) architecture is shown schematically. The absorber/detector membranes 
are suspended in a frame of silicon. Electrical and thermal connections to the frame are made over thin silicon legs. The 
backshort/termination assembly fits into keyed structures in the detector grid frame. Electrical connections are made to 
the SQUID multiplexer by wrap-around vias that are made along the grid frame, b) A prototype 8x8 detector frame is 
shown. For this prototype, the detector leads are fanned out in the plane of the detectors rather than wrapped around to 
the back of the assembly as they will be for the 32 x 40 arrays. 71 

be mounted to the SCUBA-2 32x40 two-dimensional time-domain SQUID multiplexer developed by NIST. 73 
Figure [5ji shows a schematic of the BUG architecture. 

The suspended silicon membranes on the front piece of the BUG assembly will be ion-implanted to achieve 
a surface impedance of 157 SI per square. This process has been demonstrated for the Atacama Cosmology 
Telescope (ACT) 74 detectors. The back short assembly, which for single-band implementations is made reflective, 
will be coated with an absorbing film to prevent unwanted resonances. This absorbing strategy is similar to that 
used on HAWC, and will result in a 50% efficiency across all four passbands. As in the ACT detectors, a TES 
consisting of a molybdenum-gold bilayer with noise-reducing normal metal bars will be fabricated at the edge of 
each pixel. 

NASA/Goddard has developed an 8 x 8 pixel version of this array as a step towards the realization of kilopixel 
arrays. One such array is shown in Fig. (5Jd. In this device, the TES leads are routed on the top layer of the 
device to bond pads at the edge of the structure. Multiplexing is done external to the detector chip. An 8x16 
version with 2 mm pixels has been fabricated for the GISMO (Goddard-IRAM Superconducting (2)-Millimeter 
Observer) Instrument. 75 
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